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We investigate the decoherence of histories of local densities for linear oscillators models. It 
is shown that histories of local number, momentum and energy density are approximately deco- 
herent, when coarse-grained over sufficiently large volumes. Decoherence arises directly from the 
proximity of these variables to exactly conserved quantities (which are exactly decoherent), and 
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not from environmentally-induced decoherence. We discuss the approach to local equilibrium and 
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I. INTRODUCTION 



In a large and possibly complex quantum system, which dynamical variables naturally 
become classical for a wide variety of initial states? This question belongs to the general 
issue of emergent classicality from quantum theory, and has recently received a considerable 
amount of attention (see Ref. for an overview). There are a number of different approaches 
to it, but common to most of them is the demonstration of decoherence: that certain types 
of quantum states of the system in question exhibit negligible interference, and therefore 
superpositions of them are effectively equivalent to statistical mixtures. 

Decoherence has principally been demonstrated for the situation in which there is a 
distinguished system, such as a particle, coupled to its surrounding environment 
More generally, we may expect that decoherence comes about when the variables describing 



rether or not this 
for a discussion 



the entire system of interest naturally separate into "slow" and "fast" , w 
separation corresponds to, respectively, system and environment. (See Ref. 
of the conditions under which the total Hilbert space may be written as a tensor product 
of system and environment Hilbert spaces). If the system consists of a large collection of 
interacting identical particles, as in a fluid for example, the natural set of slow variables 
are the local densities: energy, momentum, number, charge etc. These variables, in fact, 
are also the variables which provide the most complete description of the classical state 
of a fluid at a macroscopic level. The most general demonstration of emergent classicality 
therefore consists of showing that, for a large collection of interacting particles described 
microscopically by quantum theory, the local densities become effectively classical. Although 
decoherence through the system-environment mechanism is expected to play a role since the 
collection of particles are coupled to each other, it is of interest to explore the possibility 
that, at least in some regimes, decoherence could come about for a different reason. Namely, 
because the local densities are almost conserved if averaged over a sufficiently large volume 
j|. Hence, the approximate decoherence of local densities would then be due to the fact 
that they are close to a set of exactly conserved quantities, and exactly conserved quantities 
obey super select ion rules. 

We will approach the question using the decoherent histories approach to quantum theory 
DDDDQD This approach has proved particularly useful for discussing emergent 
classicality in a variety of contexts. In particular the issues outlined above are most clearly 
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expressed in the language of decoherent histories. The central object of interest is the 
decoherence functional, 

D(a,a!) = Tr (P an (t n ) ■ ■ ■ P ai {h) pP^h) ■ ■ -P<(*„)) ( l ) 

The histories are characterized by the initial state \%) and by the strings of projection 
operators P a (t) (in the Heisenberg picture) at times t\ to t n (and a denotes the string of 
alternatives a± ■ ■ ■ a n ). Intuitively, the decoherence functional is a measure of the interference 
between pairs of histories a, a'. When it is zero for a ^ c/, we say that the histories are 
decoherent and probabilities p(a) = D(a, a) obeying the usual probability sum rules may be 
assigned to them. One can then ask whether these probabilities are strongly peaked about 
trajectories obeying classical equations of motion. For the local densities, we expect that 
these equations will be hydrodynamic equations. 

We are generally concerned with a system of N particles interacting through a potential 
and are therefore described at the microscopic level by a Hamiltonian of the form 

tf = £f| + £^K-q*) (2) 

j e>j 

We are particularly interested in the number density n(x), the momentum density g(x) and 
the energy density /i(x), defined by, 

n(x) = £ *(x - q,-) (3) 

j 

g(x) = E Pi <$(x-qj) (4) 

3 

h(x) = E " q,") + E Vjtiqj - qe)6{x - q,-) (5) 

We are generally interested in the integrals of these quantities over small volumes, which will 
have the effect of smearing out the (^-functions. Integrated over an infinite volume, these 
become the total particle number N, total momentum P and total energy H, which are 
exactly conserved. It is also often useful to work with the Fourier transforms of the local 
densities, 

n(k) = E e <k * (6) 

3 

g(k) = e P, ^ (7) 

3 

m = E + E v je (^ - q*) ^ (8) 

3 zm e>j 



These quantities tend to the exactly conserved quantities in the limit k = |k| — > 0. 

There is a standard technique for deriving hydrodynamic equations for the local densities 
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. It starts with the continuity equations expressing local conservation, which have 



the form, 

§ + V-J = (9) 

where a denotes n, g or h (and the current j is a second rank tensor in the case of g). It 
is then assumed that, for a wide variety of initial states, conditions of local equilibrium are 
established after a short period of time. This means that on scales small compared to the 
overall size of the fluid, but large compared to the microscopic scale, equilibrium conditions 
are reached in each local region, characterized by a local temperature, pressure etc. which 
vary slowly in space and time. Local equilibrium is described by the density operator 
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Z~ l exp (- J d 3 x P(x) \h(x) - #(x)n(x) - u(x) • g(x)]) (10) 



where (3, p and u are Lagrange multipliers and are slowly varying functions of space and 
time. (3 is the inverse temperature, u is the average velocity field, and p is related to 
the chemical potential which in turn is related to the average number density. (Note that 
the local equilibrium state is defined in relation to a particular coarse-graining, here, the 
anticipated calculation of average values of the local densities. Hence it embraces all possible 
states that are effectively equivalent to the state Eq.([TU]) for the purposes of calculating those 
averages.) The hydrodynamic equations follow when the continuity equations are averaged 
in this state. These equations form a closed set because the local equilibrium form depends 
(in three dimensions) only on the five Lagrange multiplier fields f3, p, u and there are exactly 
five continuity equations (jHJ) for them. (More generally, it is possible to have closure up to 
a set of small terms which may be treated as a stochastic process. See Refs. , for 

example.) 

The decoherent histories approach to quantum theory offers the possibility of a much more 
general derivation of emergent classicality than that entailed in the standard derivation 
of hydrodynamics. The standard derivation is rather akin to the Ehrenfest theorem of 
elementary quantum mechanics which shows that the averages of position and momentum 



operators obey classical equations o 
must involve much more than that 



motion. Yet a description of emergent classicality 
lj. Firstly, it must demonstrate decoherence of the 



local densities, thereby allowing us to talk about probabilities for their histories. Secondly, 



it should not be restricted to a special initial state. Whilst it is certainly plausible that 
many initial states will tend to the local equilibrium state, the standard derivation does 
not obviously apply to superpositions of macroscopic states, which are exactly the states a 
description of emergent classicality is supposed to deal with. 

This paper is part of a general programme, initiated in Refs.^, 15, to obtain a more 
general derivation of hydrodynamic equations from the underlying quantum theory, using the 
decoherent histories approach. The aim in particular is to consider reasonably general classes 
of initial states and to demonstrate decoherence of the local densities, without appealing to 
environmentally-induced decoherence, and to show that the probabilities for histories are 
peaked about equations of motion of the hydrodynamic type. In this paper the particular 
system we will apply the programme to is a chain of linearly coupled oscillators. 

The general sketch of the programme, which we will work out in detail in this paper, is 
as follows. We start from the simple observation that exactly conserved quantities define 
an exactly decoherent set of histories, essentially because the projectors in the decoherence 
functional commute with the Hamiltonian ^|. It is therefore expected that the histories 
will remain approximately decoherent as we go from k = to non-zero values of k in the 
local densities n(k), g(k), h(k). In Ref.fl^ it was shown that a useful way to organize this 
idea is to decompose the initial state of the system into a superposition of states, |n, g, h), 
which are approximate eigenstates of the local densities. It is then very plausible (and 
verifiable in specific models) that such states remain approximate eigenstates of the local 
densities under time evolution, for sufficiently small k (since it is clearly exactly true in the 
limit k — > 0). Here, "sufficiently small" means that k~ x should be much greater than the 
correlation length of each local density eigenstate. The preservation in time of these states 
means that histories of them will be approximately decoherent. 

Given decoherence we may then look at the probabilities for histories. Decoherence 
also indicates that each element of the superposition of local density eigenstates may be 
treated separately. We therefore consider the probabilities for histories of local densities 
with the local density eigenstate as the initial state. For sufficiently coarse-grained histories 
the probabilities for the local densities are strongly peaked at each time about the average 
value of the local densities, averaged in the local density eigenstate |n, g, h). Since the local 
densities are sums of one-particle operators (to lowest order in interactions), this is the 
same as averaging in the one-particle reduced density operator pi constructed by tracing 
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\n, g, h)(n, g, h\ over all but one particle states. The density operator p\ is clearly not the 
same as the (one-particle version of the) local equilibrium state (fTUj) . although the two states 
are clearly very similar, since they are both very localized in the local densities. Hence, to 
complete the derivation of the hydrodynamic equations, it is necessary to show that pi tends 
to the local equilibrium state after a period of time. This is clearly extremely plausible on 
physical grounds and may be proved in explicit cases, as in this paper. (And indeed, this is 
much weaker than asserting that any initial state tends to the local equilibrium state.) 

In brief, the whole story works in particular models contingent only on constructing 
local density eigenstates and showing that they have the desired properties: that they are 
preserved in form under time evolution for sufficiently small k, and that they are effectively 
equivalent to the local equilibrium distribution after a period of time. The point of this 
paper is to show this for the linear oscillator chain. 

The detailed connection between conservation and decoherence is discussed in Section 2, 
as is the construction of approximate eigenstates of the local densities. 

In Section 3, we describe the dynamics of the linear oscillator chain. We consider two 
types of chain: the simple chain, where only neighbouring particles are coupled, and the 
bound chain, where each particle is in addition bound to an origin by a harmonic potential. 
We consider both finite and infinite chains. The most important results are the correlation 
functions, which establishes the scale on which coarse-graining is required for decoherence. 

In Section 4, as a preparation for proving decoherence of the local densities of the chain of 
oscillators, we consider a simplified set of variables, namely, the total momentum contained in 
a subsection of the chain. We show that the eigenstates of this quantity remain approximate 
eigenstates under time evolution as long as the size of the chain subsection is much greater 
than the correlation length. 

In Section 5, we consider the local densities of the chain. We prove that approximate 
eigenstates remain approximate eigenstates, for k~ l much larger than the correlation length 
of the chain. 

In Section 6, we consider the probabilities for histories. In the case of a finite simple 
chain, we show that the averages of number and momentum density obey a closed set of 
equations (although there is no evolution to local equilibrium in this case). For the infinite 
bound chain, we show that the density operator pi does indeed tend to the local equilibrium 
state and we derive the resultant hydrodynamic equations. 
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We discuss our results in Section 7. 

The idea that local densities should define a natural set of decoherent histories as a result 
of their approximate conservation was first put forward by Gell-Mann and Hartle j^. This 
idea, and the related possibilit y o f deriving hydro dynamic equations, has been developed by 
numerous authors (13, 15l lla . Il8l Il9l |20j]. This work is perhaps most closely related to 
that of Brun and Hartle |l8[, who analyze the linear oscillator chain using the decoherent 
histories approach. Their approach was rather different in that they considered coarse- 
grainings in which the centre of mass coordinates of chain subsections were specified, rather 
than the local densities considered here, and they evaluated the decoherence functional 
explicitly, rather than examine the evolution of eigenstates of the variables of interest. 



II. DECOHERENCE AND CONSERVATION 



conservation, It is well- 
The simple reason for 



We begin by describing the connection between decoherence anc 
known that exactly conserved quantities are exactly decoherent jl7]. 
this is that the projectors commute with the unitary evolution operator. The projectors 
P ak on one side of the decoherence functional (JTJ) may therefore be brought up against 
the projectors P a > k on the other side, hence the decoherence functional is diagonal. (In 
the situation considered here, in which there are three conserved quantities involved, these 
quantities must in addition commute with each other, but this is clearly the case.) 

We would like to extend this idea to approximate decoherence in the case of approximate 
conservation. It turns out that the above argument is better formulated in a different way 
for the purposes of generalization [14]. Suppose the initial state is pure and consider the 
decoherence functional, 

D(a,a!) = Tr (P an U n ^ n ■ ■ ■ P a2 U 12 P ai W^P^U^P^ ■ ■ • Ul_ 1>n ) (11) 

where U12 is the usual unitary evolution operator between times t\ and t 2 . Suppose the 
histories are projections onto some conserved quantity, Q. Now let the initial state be a 
superposition of eigenstates of Q, 

m = ^(\Qi)+\Q*)) (12) 

where 

Q\Qa) = Qa\Qa) (13) 



and a = 1, 2. Since the P^'s are projections onto Q, we have P a \Q a ) = \Qa), if ot is suitably 
chosen, otherwise we get P a \Q a ) = 0. Hence the only non-zero off-diagonal terms of the 
decoherence functional are of the form, 

D(a,g/) = ^Tr (P an U n ^ n ■ ■ ■ PMQJ (Q 2 \Ul 2 P a , ■ ■ ■ Ul_ 1>n ) (14) 

But Q is conserved, hence Ui 2 \Q a ) = \Q a ) an d 

D(a,a!) = \Tr (P an U ■ ■ ■ P a2 \Qi){Q2\P a > ■■■tf) (15) 

Proceeding in this way to the end of the chain, 

D(a,a!) = (Q 2 \Q 1 ) = (16) 

for all pairs of distinct histories a, a'. Hence decoherence comes about because neither the 
projections nor the unitary evolution disturb the states \Q a ), and hence the two orthogonal 
states \Q 2 ) are brought together at the final time and overlapped to give zero. 

Let us now suppose that we have some operator Q such that under time evolution, its 
eigenstates are mapped into approximate eigenstates. That is, we initially have Eq. (jl3|) . but 
under evolution to time t, 

Q(t)\Qa) ~ (Q(t))\Qa) (17) 

(where the average on the right-hand side is in the state \Q a ))- More precisely, this can be 
expressed as 

~WW <K 1 (18) 

where 

(AQ(t)) 2 = (Q\t)) - m)f (19) 

Eq. (|T5j) means that the state remains strongly peaked in the variable Q under time evolution. 
The states are then approximate eigenstates of the projectors at each time as long as the 
widths of the projectors are chosen to be much greater than (AQ(t)) 2 . The same argument 
goes through although this time only approximately. Approximate decoherence is therefore 
assured for sufficiently coarse-grained histories of operators Q and superpositions of initial 
states each of which have the property that they remain strongly peaked in Q under time 
evolution (as characterized by Eq.fJTHJ)). 
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A simple example is the case of the coherent states of the harmonic oscillator. These states 
are preserved in form under time evolution, hence will always be approximate eigenstates 
of projections onto position, momentum, or phase space, provided that the widths of the 
projections are chosen to be much greater than the uncertainties in p and q in the coherent 
states. In this example, there is no obvious local conservation law. For this reason, it 
is perhaps more accurate to speak of approximate determinism, rather than approximate 
conservation. So very broadly speaking, approximate decoherence of histories will arise when 
there is an approximate determinism in the underlying quantum theory. 

Returning to the local densities, we require a set of states \n, g, h) which are eigenstates 
of all three local densities. Since the local densities do not commute with each other, except 
in the limit k — > 0, we can only find states which are approximate eigenstates. The number 
and momentum density are both operators of the form, 

TV 

A = A n (20) 

n=l 

as is the local energy density, if we ignore the interaction term. For such operators it follows 
that 

(AAf = £(AA0 2 + E *( A n, Am) (21) 

and 

(A) 2 = J2(A n )(A m ) (22) 

n,m 

A state will be an approximate eigenstate of the operator A if 

W" « 1 < 23 > 

The expression for {A) 2 potentially involves N 2 terms, as does the expression for (AA) 2 , but 
the latter will involve only N terms if the correlation functions cr(A n , A m ) are very small 
or zero for n ^ m. So simple product states will be approximate eigenstates and will have 



(AA) 2 / (A) 2 of order 1/N. (See Refs. 1J, [l5| for more detailed examples this argument). 

Under time evolution, the interactions cause correlations to develop. However, the states 
will remain approximate eigenstates as long as the correlations are sufficiently small that the 
second term in Eq.(|21|) is much smaller than order iV 2 . The interactions and the subsequent 
correlations are clearly necessary in order to get interesting dynamics and in particular the 
approach to local equilibrium. The interesting questions is therefore whether there is a 



regime where the effects of interactions are small enough to permit decoherence but large 
enough to produce interesting dynamics. The fact that the variables we are interested in 
are locally conserved indicates that there is such a regime. The important point is that 
the local densities become arbitrarily close to exactly conserved quantities as k — > 0. This 
means that, at any time, (AA) 2 / (A) 2 becomes arbitrarily close to its initial value (which is 
of order 1/N) for sufficiently small k. 

In the examples we look at it in the following sections, we will see that an uncorrelated 
initial state develops correlations with a typical lengthscale (or extending to a certain number 
of particles down the chain). These correlations typically then decay with time. What we 
will find is that the second term in Eq. (J21|) will remain small as long as k~ l is much greater 
than the correlation length. Hence the key physical aspect is the locality of the interactions, 
meaning that only limited local correlations develop, together with the coarse-graining scale 
k~ l which may be chosen to be sufficiently large that the correlation scale is not seen. 
Differently put, as k increases from zero, departing from exact decoherence, it introduces 
a lengthscale k" 1 . Since the decoherence functional is a dimensionless quantity, clearly 
nothing significant can happen until k~ l becomes comparable with another lengthscale in 
the system. The natural scale is the correlation length in the local density eigenstates. 

The scheme described here would be executed most transparently if we used states which 
become exact eigenstates of the conserved quantities in the limit k — > 0, thereby always 
maintaining the closest connection with exact conservation. In the next section we will 
in fact use Gaussians as the approximate eigenstates, because they are the easiest states 
to work with. These will not be exact eigenstates of the exactly conserved quantities in 
the k — > limit, although this not in fact matter very much, for reasons outlined above. 
(Furthemore, the decoherence functional is always exactly diagonal for any initial state in 
the k — > limit, for the reasons stated at the beginning of this section, but we do not need 
to exploit this here.) 

III. CHAINS OF OSCILLATORS 

In this and the following sections, we show how the general programme outlined above 
may be worked out in detail in linear oscillator models. These have the advantage that they 
can be solved exactly. In particular, the time-development of the correlation functions and 



10 



eigenstates of the local densities can be computed reasonably explicitly. 



A. The Models and Their Classical Solutions 

We consider a chain of point particles which are coupled to each other by a nearest- 
neighbour linear interaction. We also allow the possibility that each particle is harmonically 
bound to one of a series evenly distributed points, separated by distance b, say. The Hamil- 
tonian is 



N 



n=l 



(24) 



where b n = nb. We will consider the two cases K = (the simple chain) and K ^ (the 
harmonically bound chain). In the bound chain case, it is also useful to consider the case 
b n = 0, which corresponds to the situation in which the whole chain moves in a harmonic 
potential. (In fact, for the classical solutions, b n is readily absorbed into q n , but this makes 
a difference to the local densities considered below). We initially consider a finite number 
TV of particles but we also consider the case of N infinite. 
The equations of motion are 



mq n + K(q n - b n ) = y 2 {q n +i - 2q n + q n . 



(25) 



where we take qN+i = Qi- This system has been discussed and solved in many places. A 



particularly useful reference for the case o 



and Robertson 
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221 ]. (See also Refs. 
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an infinite chain is the treatment by Huerta 



24, 



25l |26j). The solution may be found by 



introducing the normal modes, Q c 



N 2man/N 



which obey 



where 



q-n b n -\- ^ ' i Qa 



Qa + wlQ a = 



K Aw . 2 / 77 < I 



H sin 

m m 



N 



The solution may be written, 

q n (t) = b n + J2 



A ' /r-n(t)9r(0) + ^=^Pr(0) 



r=l I 



(26) 



(27) 



(28) 



(29) 
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where 



1 N 

/»(*) = - Y.z 2man ' N cos(u a t) (30) 



$ n (t) _ 1 f wjvSinM 

IT " iv^ e ~ (31) 

Here, Q 2 = (K + 2u 2 )/m. The solution for p n (t) is given by 

p n (t) = mq n (t) (32) 
In the limit of an infinite number of particles the solution is 

oo 

q n {t) = K + 



fr-n(t)q r (0) + ^p-Pr(0) 
mil 



(33) 



For the simple chain, K = 0, the solution is then given in terms of Bessel functions [21|, and 
we have 

f r (t) = J 2r (2ujt) (34) 

where uo 2 = u 2 /m, and 

g r (t) = VI [*dtf J 2r {2ujt') (35) 
Jo 

The appearance of the time integral in the expression for g r (t) is in fact related to the motion 
of the centre of mass of the whole chain, which is not of interest for our considerations. (We 
imagine the whole system is contained somehow so that there is no wavepacket spreading). 
It also somewhat obscures the discussion of correlations, which is the main thing we are 
interested in. The relevant behaviour in g r (t) is best exhibited in terms of the difference 
variables, g r+ \ — g r . Using a simple recurrence relation for the Bessel functions, these are 
given by 

g r+1 (t) - g r (t) = -2fiJ 2r+1 (2u;t) (36) 

We will discuss this in more detail below. 

For the bound chain, K ^ 0, it is most useful to work in the regime in which the 
interaction between particles is much weaker than the binding to their origins, so v 2 << K. 
In this case, the solution then is [21 1 . 

f r (t) w Jr^Vtt) cos (ttt - nr/2) (37) 

and 

g r (t) w Jri^flt) sin (fit - vrr/2) (38) 
12 



where 7 = (uj/Q) 2 , so 7 << 1. 

The general behaviour of the solutions in both cases is easily seen. The functions / r _ n (i) 
and g r - n {t) loosely represent the manner in which an initial disturbance of particle r affects 
particle n after a time t, and is given in both the bound and unbound case by the properties of 



Bessel functions. Recall the following forms of the Bessel functions 27[ • For small arguments 
we have 

J »w^(I)"(r(^-TO) + '--) < 39 > 

(This is for n > 0. For n < we use J_ n (x) = (— l) n J n (x)). For large arguments we have 
the asymptotic form 

J n (x) ~ — cos (x - vm/2 - tt/4) (40) 



\TTX, 

Hence the Bessel functions start out at zero (except for n = 0, where Jo(0) = 1), grow 
exponentially, and then go into a slowly decaying oscillation. For large n and fixed x we 
have 

j n ( X ) ~ (2™)-^ (^y (41) 

A point not immediately obvious from these standard asymptotic forms, and which will 
turn out to be important, is that the different Bessel functions each go into their oscillatory 
regions at different values of x. In particular, one can estimate from the plots of the Bessel 
functions that J n (x) goes into its oscillatory regime when x is of order n, by which stage 
J n (x) is therefore of order n -1 / 2 . In terms of the behaviour of the chain, this means that 
distant pairs of particles never come to influence each other very much, even after long 
periods of time: at short times, particles separated by n have exponentially suppressed 
correlation, like x n /n\ and at long times, their correlations are also suppressed, like rT x l 2 . 
This particular aspect turns out to be crucial for our purposes. 

Another important observation is that in the oscillatory regime, the Bessel function J n (x) 
has only a very limited dependence on n, namely it has the form (J4(Jj) for some n, plus the 
three possible phase shifts of tt/2. The significance of this for the chain is that when the 
functions f r - n (t) and g r _ n {t) have entered the oscillatory regime, the conditions at particles 
r and n and everywhere in between are approximately the same. This feature is clearly 
relevant to the approach to local equilibrium. 
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B. Correlation Functions 



As described in Section 2, we are interested in the time-development of eigenstates of the 
local densities, and this boils down to the behaviour of the various correlation functions of 
the system. We define the correlation function 

a(A,B) = ±(AB + BA)-(A)(B) (42) 

Because the system is linear the classical solutions described above may be used to discuss 
the solutions in the Heisenberg picture in the quantum case. In fact, the only quantum 
calculations in this paper have essentially been done already in Section 2, and the remaining 
analysis is essentially the same as for a classical stochastic system. 

For simplicity, we will concentrate on Gaussian initial states (which will of course remain 
Gaussian, because they system is linear), although these will be sufficient for our purposes. 
The variances of these Gaussians are restricted by the requirement that the state is pure. 
(They are also of course restricted by the uncertainty principle). We will consider two 
different types of Gaussian initial states which can be approximate eigenstates of the local 
densities. The first type we consider are product states, so have no initial correlations 
between different particles. The second type are the coherent states of the normal modes, 
and are naturally suggested by the normal mode decomposition Eq.(|26p. We will see in 
the next section that these are eigenstates of the local densities as long as the correlation 
functions remain sufficiently small. 



1. Normal Mode Coherent States 



Taking the second type first, we therefore consider the set of Gaussian states which have 

h 

2fflU n 
1 

2 

a(Q a ,K ) = (45) 



<Q^QD = t^-<W (43) 

a(K ai K* p ) = -hmwJap (44) 



where K a is the momentum conjugate to Q a and note that Q* a = Q- a . Because these are 
the coherent states of the harmonic oscillator, these correlation functions all remain of this 
form under time evolution, and the only time-development of the states is in terms of their 
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centres, Q a (t), K a (t), which follow the classical equations of motion. For the case of the 
simple chain, K = 0, we have u a = when a = N. This corresponds to the centre of mass 
of the whole chain, and may be quite simply omitted. 

In terms of the original coordinates q n and p n , we have the correlation functions 

•(-•J^Si^ (46) 

i N 1 

°(Pn,Pm) = ^ E 2 kmUJa ^^'^ (47) 
a=l 

<?(q n ,Pm) = o (48) 

These correlation functions are constant in time, and this feature makes this case a useful 
one to study. Also notice that (Aq n ) 2 and (Ap n ) 2 are independent of n. The correlation 
functions will typically decay very rapidly with increasing \n — m\, since they are sums of 
rapidly oscillating terms. This is especially true in the case K ^ with K » u 2 , since 
then u a , Eq. (j28j) . is a constant to leading order. (We will a similarly effect in more detail in 
the next subsection). In the case K — 0, we have 

"" = S sm (lvJ (49) 

The correlation function a(p n ,p m ) clearly still decays for large \n — m\, but this is less 
obvious for a(q n , q m ), Eq. (J46|) . because the denominator becomes very small close to a = N. 
However, if we eliminate by hand a small cluster of modes close to a = N, we get satisfactory 
decay properties for this correlation function, and we will assume that this has been done. 

2. Uncorrelated Initial States for the Infinite Chain 

We now consider uncorrelated initial states for the infinite chain, so the correlation func- 
tions o~(q n , q m ), cr(q n ,p m ) and a(p n ,p m ) all vanish at the initial time for n ^ m. We are then 
interested in the behaviour of these three types of correlation functions at later times. 

From the solution Eq.()33J). with the assumption of no initial correlation between the 
particles, we have 

0-(q n (t),q m (t)) = (/r-«W/r-m(f)(A? r ) 2 + ^9r-«W9r-m(t)(Aj) r ) 2 

+ mQ ifr-^dr-mit) + f r -m(t)g r -n{t)) a(q r ,p r )J (50) 
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and similarly for a(q n (t),p m (t)) and a(p n (t),p m (t)). These expressions simplify if we make 
the further assumption that the initial values of (Aq r ) 2 , (Ap r ) 2 and a(q r ,p r ) are independent 
of r (we will show below how to go beyond this assumption). In this case we obtain 

v{qn{t),q m {t)) = a nm {t){Aq) 2 + 2e nm (t)cr(q,p) + d nm {t){Ap) 2 (51) 
(j(q n (t),p m (t)) = b nm (t)(Aq) 2 + [a nm (t) + k nm (t)]a(q,p) + e nm {t){Ap) 2 (52) 
a(p n (t),p m (t)) = c nm (t)(Aq) 2 + 2b nm (t)a(q,p) + a nm (t)(Ap) 2 (53) 

where 

a n m{t) = J2fr-n{t)fr-m{t) (54) 
r 

b nm (t) = mJ2fr-n(t)f r -rn(t) (55) 

r 

Cnm(t) = m 2 ^ fr-n(t) f r -m(t) (56) 

r 

dnmit) = } J2g r -n(t)g r - m (t) (57) 

m z \l z r 

e nm (t) = — 7yj^2 9r-n(t)9r-m(t) (58) 
knmit) = — ^ fr-n{t)g r -n{t) (59) 

Since the coefficients f n (t) and g n {t) are all given by Bessel functions, these expressions can 
be evaluated using the following Bessel function addition theorem: 

oo 

J n (2x)= Jn-k(x)J k (x) (60) 

k=— oo 



2Qt (n — m)7T 

2 V ; 



(63) 



For the bound chain, K ^ 0, the coefficients are 

2a nm (t) = 5 nm + J n _ m (2jQt) cos 2Qt - ~(n - m)n (61) 

2b nm (t) = — mQJ n ^ m ( y 2 r yQt) sin 2Qt — -(n — m)ir (62) 

2c nm (t) = (mtt) 2 5 nm - (mVt) 2 J n _ m (2^Vtt) cos 

dnmit) = (mQ)~ 4 c nm (t) (64) 

e„ m (t) = -(mVty 2 b nm (t) (65) 

fcnmW = -{mVty 2 c nm (t) (66) 

All of these coefficients, and hence all of the correlation functions, decay exponentially for 
large \n — m\. Furthermore, in the limit t — > oo we have 

<j{q n (t), q m {t)) -> \b nm {{Aqf + -I^(Ap) 2 ) (67) 
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a(q n (t),p m (t)) -> (68) 
<r(p n (t),p m (t)) -> i(5 nm (m 2 fi 2 (Ag) 2 + (Ap) 2 ) (69) 

We also find that (q n (t)) — > b n and (p n (t)) — * as t — > oo. We thus obtain an equilibrium 
distribution, since the corresponding phase space distribution function for the whole chain 
is 



n 



kT 

where we identify the temperature as 



PL + - m n 2 (q n -b n ) 2 
2m 2 W ; 



(70) 



fcT = — f m 2 n 2 (Aq) 2 + (Ap) 2 ) (71) 
2m v 7 

For the simple chain, K = 0, we have 21], 

2a nm (t) = 5 nm + J 2 n-2m(4wt) (72) 

2b nm (t) = mUJ [J 2 n-2m-l(4^t) - J 2 n-2m+i(4u;t)] (73) 

2c nm (t) = (mw) 2 [J 2n - 2m+2 (4^t) + J 2 n-2m-2 (4w*) - 2J 2n _ 2m (4^t) (74) 

+ 25 nm ] (75) 



Mart I 1 Mwt 

/ o ^ J^fe)*, (76) 

= 1 / J2n-2m(y)dy (77) 



2a;m : 



4u;m jo 

(The explicit form of k nm (t) is not required). The first three of these coefficients, like the 
bound case, are exponentially suppressed for large \n — m\. This means that the behaviour 
of the correlation cr(p n (t),p m (t)), which depends only on these three coefficients, has the 
expected behaviour, but the other correlation functions do not have this property. In par- 
ticular, from the behaviour of cr(q n (t),q m (t)), the whole chain appears to become highly 
correlated. In the long-time limit, we find 

a(q n (t),q m (t)) -> ^6 nm (Aq) 2 + ^-^(Ap) 2 (78) 

cr(q n (t),p m (t)) - ^( A ^) 2 ( 79 ) 

1 1 
a(p n (t),Pm{t)) -> ^(mu) 2 [2S nm - 5„, m _i - 5„, m +i] (Ag) 2 + -6 nm (Ap) 2 (80) 

(where for simplicity we have taken cr(q n ,p n ) = initially). This is not an equilibrium 
distribution, and the behaviour of (Aq n (t)) 2 at late times indicates diffusive behaviour. 
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This growth of correlations and variances is essentially unphysical, since any realistic 
system is contained in some way, so the spreading cannot proceed beyond the size of the 



container (and indeed there is no such effect in the if / case). It is, however, difficu 



to factor out this unphysical effect in a convenient way. The long-time limits of Ref. 



2l|, 



t 



which we have followed closely, are dominated by this spreading effect in the K = case. In 



an attempt to understand this, the authors of Ref. [2l| considered a different set of models 



in Ref. |22|. in which the chain was fixed at one end. This avoided the diffusive growth 
encountered above, but still led to significant correlations between all particles on the chain 
in the long-time limit. This in turn spoils the desired behaviour of the local densities, 
discussed below. The upshot of this is that it is not possible to prove decoherence of the 
local densities in the case of the infinite chain with K = 0. 

There are two simple ways in which the above results are easily generalized. First of all, 
note that although we are focusing on Gaussian initial states, expressions for the correlation 
functions such as Eq. (|50|) are in fact valid for any initial state, because of the linearity of 
the dynamics. 

Secondly, in deriving Eqs. (j51)) - (j5H|) . we assumed that the initial variances are independent 
of r. This assumption was necessary in order to be able to apply the Bessel function 
addition theorem Eq.()60p. and thereby obtain explicit expressions for the coefficients, a nm , 
b nm etc. This is too restrictive, since it means that certain hydrodynamic variables (such as 
temperature, which depends on (Ap r ) 2 ), are obliged to be constant along the chain. One 
can see, however, that these results easily extend to the case in which the initial variances 
(Ag r ) 2 , (Ap r ) 2 and a(q r ,p r ) have a slow dependence on r along the chain. The point is that 
because of the decay of the functions f n - r , g n -r for large \n — r\, the sum over r in Eq.(jS0J) is 
effectively restricted to a finite range, namely, over the (actually quite small) range for which 
significant correlations exist. As long as the initial variances vary significantly with r only 
on a range larger than the correlation range, then the calculation of correlation functions 
is effectively equivalent to the case in which the variances are completely independent of r. 
This means that in place of Eq. (j51j) . for example, we get the more general result 

cr(qn(t),q m {t)) ~ a nm (t)(Aq r ) 2 + 2e nm (t)a(q r ,p r ) + d nm (t) (Ap r ) 2 (81) 

where r on the right-hand side is taken to be mid-way between n and m, and for n = m we 
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get 

(Aq n {t)) 2 « a nn {t){Aq n ) 2 + 2e nn {t)a{q n ,p n ) + d nn {t){Ap n ) 2 (82) 

This simple observation is important for obtaining interesting hydrodynamic equations be- 
cause it allows for the possibility of the system tending towards local equilibrium, rather 
than equilibrium of the whole chain. 



IV. COARSE GRAINING BY CHAIN SUBSECTIONS 

Although we are ultimately interested the local densities for the chain variables, we will 
first consider some simpler variables whose analysis is highly instructive. Namely, we take 
the variables of interest to be the total momentum in a subsection of the chain containing 
M particles, so we define 

M 

P M = Y,Pn (83) 
n=l 

where M <C N. This is not quite the same as a locally conserved quantity, but it is 
very similar, since, for the simple chain, the total momentum is conserved. Therefore Pm 
is an exactly conserved quantity when M = N, and otherwise we might expect it to be 
approximately conserved. 

As outlined in Section 2, to show that these variables are approximately decoherent, we 
need to show that initial states satisfying the condition 

« i «*) 

will continue to satisfy it under time evolution, hence the initial state remains an approximate 
eigenstate. We have 

M M 

(AP M ) 2 = £ E °(Pn,Pm) (85) 
n=l m=l 

and 



M 

(PM) 2 = Y,(Pn)(Pm) (86) 

n,m 

We take the case in which the initial correlations are zero and the initial variances are the 
same all along the chain, and we also take the initial value of cr(q,p) to be zero. We thus 
obtain 

(AP M (t)) 2 = C M (Aq) 2 + A M (Ap) 2 (87) 
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where 

M M 

A M (t) = E E a nm(t) (88) 
n=l m=l 

M M 

C M (t) = E E Cnm(t) (89) 
n=l m=l 

and the coefficients a nm (t) and c nm {t) are given by Eqs.(|72 |l . (j73j) . The two terms on the 
right are very similar, so for simplicity will will concentrate on the second one. (Note also 
that these terms do not suffer from the spreading problem discussed in the previous section. 
This is an advantage of working with momenta, rather than positions). From the expression 
(J?2J), we have 

I M M 

A M(t) -EE J2n-2m{^t) (90) 
n=l m=l 

(where we have assumed that N is sufficiently large that it is effectively equivalent to the 
iV — > oo case). Our aim is now to show that Am <§C M 2 for all times, for then the condition 
Eq.flSID will be satisfied. 

The expression for Am cannot be evaluated exactly, but its properties may be seen from 
the asymptotic forms of the Bessel functions (|39 j) . ([40 )1 . For small times, when the Bessel 
functions are all in the exponential regime, J2n-2m is exponentially suppressed for large 
| n — m\. For larger times, the Bessel functions start to go into their oscillatory form, where 
they are already small. Furthermore, because in the oscillatory regime they depend on 
n — m only through a simple phase, most of the terms in the sum over n and m cancel 
out. Proceeding along these lines one can see that Am will not come anywhere close to M 2 
except for small values of M. These features are easily confirmed by plotting Aj^if) for 
different values of M. For example, with M — 5, quickly decays to a value of about 

0.1, clearly much less than M 2 = 25, thereafter going into a slowly decaying oscillation. 

As we shall see, the variances of the local densities differ from the momentum of a 
chain subsection in that they are more complicated functions of the correlation function, for 
Gaussian states, but the physical understanding of their behaviour is essentially the same, 
which is why this simple example is instructive. 
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V. DECOHERENCE OF LOCAL DENSITIES 



We now come to the main point of this paper, which is to examiner the behaviour of the 
local densities for the oscillator chain. They are 



N 



n\x) 



9{x) 



h(x) 



n=l 

N 

^2Vn^{q n - x) 



71=1 

N 

E 

71=1 



+ y (9n - Qn-i) 2 + \K{q n - b n ) 2 



S{q n - x) 



(91) 
(92) 
(93) 



Again it will often be very useful to work with the Fourier transformed local densities, 

N 

n(k) = $> ifc9n (94) 

n=l 
N 

9(k) = Ep/ 9 " (95) 

71=1 

N In 2 v 2 1 1 

Kk) = E ir L + ^(ln-q n -i) 2 + -K(q n -b n ) 2 e ik «« (96) 



71=1 L 



The local number and local energy density are locally conserved. The local momentum 
density is locally conserved only for the case of the simple chain, K = 0. 

We consider states which are approximate eigenstates of the local densities. We will use 
Gaussian initial states, and we expect that these will be approximate eigenstates of the local 
densities if we choose the correlation functions a(q n , q m ), cr(q n ,p m ) and <j(p n ,Pm) to be zero, 
or at least sufficiently small, for n ^ m. 

For the Gaussian states we consider here, the computation of variances of the local 
densities is facilitated by the identity, 



exp [a n (q n - (q n )) + f3 n (p n - (p n ))] J) 



cxp 



E 

nj 



n,Pj , 



(97) 
(98) 



All of the variances of interest will therefore be functions of the three basic types of corre- 
lation functions cr(g n ,g m ), a{q ni p m ) and cr(p n ,p m ) discussed in Section 3, and the physical 
discussion will in fact be very closely related to that of the simple case discussed in Section 
4. Actually, the formula Eq. (j9*H|) also holds to quadratic order in a n and (3 n for any state, 
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so the results derived below on the basis of Gaussian states will be valid for arbitrary states 
for small k. 

We consider first the Fourier transformed number density n(k). In a general Gaussian 
state, we have 

(n(k)) = J>**> = £exp ( ik{qj) _ k 2 (A qj f) (99) 

3=1 3=1 \ ^ / 

and 



(An(fc)) 2 = (n\k)n(k)) - \(n(k))\' 



N N 



= J2J2( eikqj )( e ~ ik9n ) (e k2 ° {q} ' qn) -l) (100) 

jr = l n=l 

The latter is to be compared with 

N N 

IWfc))l 2 = EE(^)(^") (ioi) 

3=1 n=l 

With an initially uncorrelated state we have cr(qj, q n ) = for j ^ n and we see that 

(An(k)) 2 = J2\( eikqj )\ 2 (e kHAq])2 - l) (102) 

3 

From this we expect that 

(An(k)) 2 <^\(n(k))\ 2 (103) 

as long as /c -1 does not probe on scales that are too short (compared to Ag^), so the state 
is an approximate eigenstate. 

When there are correlations, as will arise over time, we expect that the state will still be 
an approximate eigenstate on lengthscales which are much greater than the lengthscale 
of correlation. As k increases from zero we have 

(Mfe)) 2 „ k\AXf 

\(n{k))\ 2 N 2 y 1 

where X = J2j 1j (the centre of mass coordinate of the whole chain). This will be very small 
as long as is much larger than the typical lengthscale of a single particle. (An(k)) 2 starts 
to grow very rapidly with k, and Eq. (jl(J3j) is no longer valid, when k~ l becomes less than the 
correlation length indicated by cr(qj, q n )- Hence the state is strongly peaked about the mean 
as long as the coarse graining lengthscale k~ l remains much greater than the correlation 
length of the local density eigenstates. 
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Consider now the local momentum density. We have 

(9(k)) = £((?;> + ik(T(q j ,p j ))exp (ikfa) - ;U 2 (A ?i ) 2 ) (105) 

3 

and at some length, we find 

(Ag(k)) 2 = J2(e ikqj ){e- ikqn ) (A jn + B jn + C jn ) (106) 

jn 

where 

Ajn = vti>j,Pn)e k2a{qi « n) (107) 
Bjn = «p„) - ik<r(q n - q v Pn)) ((Pj) + ika( qj - q n , Pj )) (e fcV <*'**> - l) (108) 
C jn = -ik ((p n ) - ika(q n ,p n )) - ik ((pj) + ika(q j ,p j )) + k 2 a(q j ,p n )a(q n ,p :j ) (109) 

We require (Ag(k)) 2 to be small in comparison to 

|<0(*)>| 2 = $>^><e-^> ({pj)(p n ) + ik(p n )a(qj,pj) 

jn 

-ik{p j )a(q ni p n ) + k 2 a(q j ,p j )a(q n ,p n )^ (110) 

Despite the complexity of these terms the interpretation is reasonably simple. As k — ► 0, 

(A<?(fc)) 2 (AP) 2 

where P = J^j is the total momentum, and it is easy to confirm that this is small (typically 
order 1/N) for the states we are using. As k increases from zero, (Ag(k)) 2 will grow, and will 
potentially contain of order N 2 terms, the same as \ (g(k))\ 2 . By inspecting the three terms 
Eqs. fll07j) - (jll0j) . however, one can see that each of them are prevented from generating order 
A" 2 terms as long as, respectively, the correlation functions cr(pj,p n ), a(qj, q n ) and u(qj,p n ) 
are suppressed for j ^ n. As in the case of the number density, (Ag(k)) 2 /\(g(k))\ 2 will start 
to grow appreciably as k~ l approaches the lengthscale indicated by cr(qj,q n )- 

Consider now the energy density. The computation of (Ah(k))) 2 is rather complicated, 
but since we are considering Gaussian states, the final conditions on the correlation functions 
will be very similar to those on the variance of the momentum density considered above, 
so we will not carry out the computation explicitly. Instead, we consider a simpler special 
case. Take the case where each oscillator is fixed to an origin at bj, and suppose that the 
binding to it is so strong that each particle is well- localized around bj. Then the integral 
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over a volume of the energy density h(x) is then approximately equivalent to taking a coarse- 
graining consisting of the energy contained in a chain subsection, similar to Section 4. We 
therefore consider the variable 

M 

h M = J2 h J ( 112 ) 



where 



(neglecting the interaction term). Following the general discussion of Section 2, we will have 
(Ah M ) 2 <C (h M ) 2 provided that 

J2^h 3 ,h n ) < (h M ) 2 (114) 

jn 

so the left-hand side must be much smaller than order M 2 . The correlation function a(hj, h n ) 
is constructed from terms like cr(<?|, q n ) and similar functions, and we have 

a(q 2 ,q 2 l )=2a(q J ,q n ) 2 (115) 

(recalling that (qj) ~ bj since the particles are tightly bound). The discussion is then very 
similar to the case of Section 4, with essentially the same result, which is that the state is 
an approximate eigenstate as long as M » 1. 

We have now shown that Gaussian states will be approximate eigenstates of the local 
densities for k sufficiently small compared to the correlation length determine by cr(qj,q n ), 
and provided also that the other correlation functions cr(qj,p n ) and cr{pj,p n ) are small for 
j 7^ n. (These results also hold for non-Gaussian states for small k). The results of Section 2 
show that all the correlation functions (excluding the infinite chain in the K — case) have 
the desired properties for all time. This proves the desired result that initial eigenstates 
of the local densities remain approximate eigenstates under time evolution, for sufficiently 
small k. 



VI. HYDRODYNAMIC EQUATIONS FOR THE LOCAL DENSITIES 

We have shown in a variety of circumstances that eigenstates of local densities are approx- 
imately preserved in form under time evolution, on sufficiently coarse-grained scales, and 
therefore superpositions of them define decoherent sets of histories. We may now therefore 
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look at the diagonal elements of the decoherence functional, representing probabilities for 
histories of these variables. These probabilities are peaked about the average values of the 
local densities, averaged in the approximate eigenstates we have been considering. (This is 
reasonably obvious, but a more detailed proof appears in the appendix of Ref.[l5j].) We will 
now show that these average values obey hydrodynamic equations. 

A. Local Conservation Equations 

The local densities satisfy the local conservation laws, 

ftW = -si (116) 

Or 

g(x) = -—-Kxn{x)+KY,b j S{q j ~x) (117) 

fc(x) = -| (118) 

They are actually more usefully written in momentum space, 

n(k) = %(k) (119) 
g(k) = tkr(k) - Kj2(q 3 - bj)e ik ^ (120) 

3 

h(k) = ikj(k) (121) 
where the currents r(k) and j(k) are given by 

r(k) = £^+* 2 £^ jr J ~ (122) 

3 m j %K 

m = £ | (f| + \K{<13 - hf + \»\<l3 - *-i) 2 ) e*« (123) 



J 



p. (e ikq i - e ikq J +1 ^ 



+ ^e>+.- ^ « (124) 

(These are clearly finite as k — > 0.) These currents are rather complicated in configuration 
space, except in the case where we neglect the interaction term, when they are given by 

2 

r(x) = £ tLSfa - x) (125) 

fa) = £ | (f| + ^te - & .) 2 ) % - *) ( 126 ) 
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Equations (|116|) - ([118p do not in general form a closed system, so do not lead to hydro- 
dynamic equations. To get a closed set, it is necessary to average these equations in a set of 
states depending on just three fields, thereby obtaining three equations for three unknowns. 
In the standard approach to deriving hydrodynamics, the local equilibrium state is chosen. 
We will discuss this below in Section C, but first we consider the simpler and instructive 
case of the normal mode coherent states. 



B. Hydrodynamic Equations in the Case of Normal Mode Coherent States 

We consider the simple chain, K = 0, in the normal mode coherent states. As we 
have shown, these states are strongly peaked in the local densities so such states define a 
decoherent set of histories. The correlations functions of these states are constant in time, 
so we do not expect a settling down to local equilibrium. However, it turns out that the 
averages of the local densities still obey a simple set of hydrodynamic equations, and this 
case turns out to be particularly transparent and instructive. Because the only dynamics in 
this case is contained in the motion of the centres (qj), (Pj), we need consider only the local 
number and momentum densities, not the energy density (which may in fact be calculated 
from the number and momentum densities in this case). Closure of the averaged conservation 
equations (JHHJ), ()120j) is obtained in this case because there are two equations and the states 
depend on just two sets of quantities, the (qj) and the (pj). 

In a general Gaussian state we have 

(n(A0) = E ex P (ikfa) - ^(A^) 2 ) (127) 
(9(k)) = E [<Pi> + ikvfapj)] exp (ik( qj ) - itfCAfc) 2 ) (128) 



j 



(r(k)) = E -((Pi) + ika( qj , Pj )f + C 



m 



2 / A „ \2 



exp [ik(qj) - -k (A^ 



3 

< 2 ' 1 



where 



2 

+ ^ E(fe+i> " 2fo> + exp (ik( qj ) - ^ 2 (A ?j ) 2 ) (129) 



Cj = —(Apj) 2 + v 2 [a(q j+1 , Qj ) - 2(Ag,) 2 + a( Qj , q^)} (130) 



Now recall that the normal mode coherent states have the following special simplifying 
features: a(qj,Pj) = 0, and the variances (Ag^) 2 , (A^-) 2 and correlation functions of the 
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form a(qj + i, qj) are independent of j. Also, we find that Cj = (which follows from taking 
the time derivative of a(qj,pj) = 0). We therefore have 

(*)) = £exp(*%,)-^ 2 (Ag) 2 ) (131) 



j 

1 



(g(k)) = exp (ik( qj ) - -k\Aq) 2 j (132) 

(r(£;)> = r p {k) + r q {k) (133) 



where 



'pW = E ^^) 2ex P - ^ 2 ( A <?) 2 ) ( 134 ) 

^(*0 = ^ E«Sf+0 - 2 fo> + ex P ( ik (*) - ^ 2 ( A ?) 2 ) ( 135 ) 

Generally, we do not expect to derive interesting hydrodynamic equations except in the long 
wavelength regime. Clearly in this case, this means k~ 2 » (Ag) 2 . Going to this regime 
(whose significance will become apparent below), and reverting back to configuration space, 
we find 

(«(*)> = E ( 136 ) 

3 

(g(x)) = EfoWfo)-*) ( 137 ) 

3 

T p ( X ) = E-^)^(fe)-^) ( 13 §) 
3 m 

The quantity r q (x) is more complicated and will be dealt with below. In these expressions, 
since we are in the long-wavelength regime, the 5-functions are to be thought of as smeared 
over a volume greater than Aq. 

It is now very useful to introduce the velocity field v(x), defined via 



(pj) = mv((qj}) = J dy mv(y)5((q j ) - y) (139) 
which inserted in Eq. (|137|) yields 

(g(x)) — mv(x)(n(x)) (140) 

This is in fact the standard relation between the velocity field and the momentum density 
11]). Inserted also in Eq. (jl38j) we obtain 

t p (x) = mv 2 (x)(n(x)) (141) 
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It is for deriving these last two equations that the long wavelength assumption is necessary. 
Most importantly, r p (x) is expressed in terms of the two fields v(x), (n(x)), which is crucial 
for closure of the equations. 

We need now to obtain a similarly simple expression for r q (x)). As it stands, Eq. (J135j) will 
not lead to a simple expression in terms of v (x) and (n(x)). To proceed further with this term 
we need to make simplifications. We are ultimately interested in deriving a wave equation for 
the number density (which one might expect on the basis of the classical equations of motion 
for the qj). The key to this is to consider small displacements of the qj about uniformly 
distributed initial values, and then to consider the fluctuation in number density about the 
constant background. We therefore write 

(q j )=jd + 5q J (142) 

where d is a constant representing the spacing between each particle and 5qj is a small 
displacement. The average number density for k~ 2 ^> (Ag) 2 and to linear order in 5qj is 
then 

(n(k)) = ^e ifc <*> (143) 

3 

= no(k) + m(A;) + • • • (144) 

where 

no(k) = X> ifei " (145) 
j 

ni(fc) = ikJ2^Qj e ikjd (146) 

and note that no(k) is constant. Inserted in the expression for r q (k), Eq. ()135|) . and assuming 
also that k" 1 ^> d, we find 

Tq (k) = dV ni (k) (147) 
This now means that in configuration space we have 

(t(x)) = mv 2 (x)(n(x)) + d 2 v 2 n x {x) (148) 

Inserting all of these results in the local conservation equations (|116)) . (jll7|) . we obtain a 
closed set of equations for v(x) and (n(x)) = uq{x) + ni(x) where 

n Q (x) =Y i 6(x-jd) (149) 

3 
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is a fixed background field (again interpreted as coarse-grained over a length of order Aq). 
Explicitly, we have 



dt 



Ft {n) = ~lL {{n)v) 

d ( d 2 v 2 s 

((n)v) = -— (n)v 2 H n x 

ox \ m 



(150) 
(151) 



The interesting special case is that in which the velocity field is small, in which case we may 
neglect the v 2 term in Eq. (jl5H) . Then combining the two equations yields the wave equation 
for nAx), 



d 2 ni 2 d 2 ni 

— c 

dt 2 dx 2 



(152) 



where c = d v jm. 



C. Local Equilibrium and Hydrodynamic Equations for the K / Case 



We now consider the case of the infinite chain with K ^ 0. In this case, we expect that 
the local density eigenstates will settle down to a local equilibrium state after a period of 
time. We will justify this important step below, but first we consider the consequences of 
a local equilibrium state, the standard assumption in derivations of hydrodynamics. This 
state is characterized by the one-particle distribution function (Wigner function) 

I fa - mv(q j )) 2 \ 



(153) 



v 2mkT(qj) J 

where /, v and T are slowly varying functions of space and time. (This is the one-particle 
distribution function for particle j - it is labelled by j since the particles are distinguishable). 
If we now average the system Eqs. (|llfij) - (|118j) . together with the currents t(x), j(x) in the 
local equilibrium state, we obtain a closed system, since we get three equations for three 
unknowns. In the case of negligible interactions and bj — 0, we find 



(n(x] 

(h(x] 
(t(x] 



Nf(x) 

mv(x)Nf(x) 

^mv 2 + hzT + \Kx 2 ^j Nf(x) 
(mv 2 + kT) Nf(x) 
^vkT + ^mv 3 ^j Nf(x) + ^ 2 (g(x)) 



(154) 
(155) 
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Inserted in Eqs. ([116)) - ()118|h the above relations give a closed set of equations for the three 
variables /, v and T. After some rearrangement, these equations are 

— + v— = -f— (159) 
dt dx dx 

dv dv 1 09 9 df Kx 

at ox m ox mj ox m 

— v— 29 9V (161) 
dt dx dx 

where 9 = kT. These are the equations for a one-dimensional fluid moving in a harmonic 
potential Note that non-trivial equations are obtained even though we have neglected 
the interaction terms in deriving them. The role of interactions is to ensure the approach 
to local equilibrium, as we discuss below. 

In these expressions, the definition of the velocity field is equivalent to Eq. (jl39j) and, 
similarly, the definition of the temperature fields is essentially equivalent to, 

£ T^^Pjf <% - x) = \kT(x)n(x) (162) 

3 

(recalling that we are working at long wavelengths, so the 5-function is coarse-grained over 
a scale of order k~ l ). Hence temperature arises not from an environment, but from the 
momentum fluctuations averaged over a coarse-graining volume. 



D. The Approach to Local Equilibrium 

Now the key point is that the states we are actually interested in are the approximate 
eigenstates of the local densities, evolved in time, or more precisely, the one-particle distri- 
bution function W\ constructed from those states. (Since W\ is the quantity that will enter 
in the computation of any averages of sums of one-particle quantities, such as the local 
densities). This is not necessarily the same as the local equilibrium distribution, Eq. (jl53j) . 
although they are clearly very similar. The averages of the local densities in the approxi- 
mate eigenstates will therefore obey the hydrodynamic equations derived above as long as 
we can show that the one-particle distribution function of these states W\ approaches the 
local equilibrium form Eq. (jl53j) after some time. 

The local equilibrium form Eq. (|153jl has a(qj,pj) = and all the other averages {pj}, (qj), 
(Agj) 2 and (Apj) 2 vary slowly in time and space (i.e. slowly with j). (Clearly a(qj,pj) has 
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to be zero or small for local equilibrium since it is the time derivative of (Aqj) 2 and (Apj) 2 ). 
Compare this with the approximate eigenstates for the case K ^ and the infinite chain. 
They are Gaussians, so their one-particle distribution function w\ is entirely determined 
by (pj), (qj), (Aqj) 2 , (Apj) 2 and a(qj,pj). From Eqs. (joT|) - (jo^|) . we see that cr(qj,pj) grows 
initially from zero, but then becomes small at late times, whilst {Aqj) 2 and (Apj) 2 approach 
a constant plus a slowly decaying factor. Moreover, (Aqj) 2 and (Apj) 2 depend on j only 
through their initial values, which as discussed (see Eq. (j82|) for example), vary slowly. In 
addition, the centres {qj(t}) and {pj(t}) go into a phase of slow time dependence and limited 
dependence on j for times sufficiently long for the Bessel functions to go into their oscillatory 
phases. These asymptotic forms are approached on a timescale ( r yil)~ 1 . Therefore, in this 
case, the local equilibrium form is indeed achieved at late times, and the average values of 
the local densities obey hydrodynamic equations. 

The final picture we have is as follows. We can imagine an initial state for the system 
which contains superpositions of macroscopically very distinct states. Decoherence of his- 
tories indicates that these states may be treated separately and we thus obtain a set of 
trajectories which may be regarded as exclusive alternatives each occurring with some prob- 
ability. Those probabilities are peaked about the average values of the local densities. We 
have argued that each local density eigenstate may then tend to local equilibrium, and a set 
of hydrodynamic equations for the average values of the local densities then follow. We thus 
obtain a statistical ensemble of trajectories, each of which obeys hydrodynamic equations. 
These equations could be very different from one trajectory to the next, having, for example, 
significantly different values of temperature. (In the most general case they could even be 
in different phases, for example one a gas, one a liquid.) 

Decoherence requires the coarse-graining scale k~ x to be much greater than the correlation 
length of the local density eigenstates, and the derivation of the hydrodynamic equations 
requires k~ 2 3> (Aq) 2 . In brief, the emergence of the classical domain occurs on lengthscales 
much greater than any of the scales set by the microscopic dynamics. 

VII. SUMMARY AND DISCUSSION 

We have shown that for the linear oscillator chain the local densities define a decoherent 
set of histories of sufficiently coarse-grained scales. The key idea is to split the initial state 
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into local density eigenstates and show that they are preserved in form under time evolution. 
The subsequent probabilities for histories are peaked about the average values of the local 
densities, and the equations of motion for them form a closed set of hydrodynamic form on 
sufficiently large scales, provided, in general, that sufficient time has elapsed for the local 
density eigenstates to settle down to local equilibrium. 

It is perhaps of interest to comment on the timescales involved. Decoherence through 
interaction with an environment involves a timescale, which is typically exceptionally short. 
Here, however, there is no timescale associated with decoherence by approximate conser- 
vation. The eigenstates of the local densities remain approximate eigenstates for all time. 
There is, however, a timescale involved in obtaining the hydrodynamic equations, namely, 
the time required for a local density eigenstate to approach local equilibrium. In this model, 
this timescale is of order (for the infinite chain in the K 7^ case). 

Another interesting general issue is the question of the relative roles of the decoherence 
through approximate conservation considered here versus decoherence through interaction 
with an environment. The point is that it is a question of lengthscales. We have demon- 
strated decoherence of the local densities starting with exact conservation at the largest 
lengthscales and then moving inwards. In this way we were able to prove decoherence with- 
out using an environment, for certain sets of histories at very coarse-grained scales whose 
probabilities are peaked about classical paths. However, in general we would like to be able 
to assign probabilities to non-classical trajectories. For example, what is the probability that 
a system will follow an approximately classical trajectory at a series of times, but then at one 
particular time undergoes a very large fluctuation away from the classical trajectory? The 
approach adopted here indicates that the probability for this is approximately zero, to the 
level of approximation used. Yet this is a valid question that we can test experimentally. It 
is at this stage that an environment becomes necessary to obtain decoherence, and indeed it 
is frequently seen in particular models that when there is decoherence of histories due to an 
environment, decoherence is obtained for a very wide variety of histories, not just histories 
close to classical. It is essentially a question of information. Decoherence of histories means 
that information about the histories of the system is stored somewhere 0, Classical 
histories need considerably less information to specify than non-classical ones, and indeed 
specification of the three local densities at any time is sufficient to specify their entire clas- 
sical histories. This is not enough for non-classical histories, so an environment is required 
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to store the information. 

Note also that "environment" need not necessarily refer to an external environment. It 
could also include the internal coordinates not fixed by the coarse-graining. These did not 
play a role in the case considered here, but would become important at finer-grained levels, 
producing fluctuations about the evolution described by the hydrodynamic equations, hand 
in hand with decoherence. This has been considered in Refs.|4l. Il3] 

Given the need for an environment at finer-grained scales, it is then of interest to ask 
whether the local densities continue to have an important role for many-body systems when 
an environment becomes necessary for decoherence. Gratifyingly, the answer is that the 
local densities, and particularly the number density, remain the naturally preferred variables 
for a many-body system coupled to an external environment, as was recently shown [29]. 
It is normally claimed that position is the preferred variable in environmentally-induced 
decoherence, but this is for a single particle coupled to an environment and is in any case 
an approximation. For a many-body system coupled to a scattering environment, with both 
described by many-body field theory, it was shown in Ref . that number density is the 
naturally decohering variable (with momentum density, as its time derivative, also entering 
in a natural way). Hence there is a smooth match with environmentally-induced decoherence 
models as we go to finer scales. 

It would be of interest to generalize to an oscillator chain with non-linear interactions. 
This is because in the linear chain, the energy in each mode is conserved, so there is no 
possibility of exchanging energy between modes, and the approach to local equilibrium is 
rather artificial. 

It would also be of particular interest to look at a gas. Many-body field theory may be 
the appropriate medium in which to investigate this, following the lead of Ref.j29|. The 
decoherent histories analysis might confer some interesting advantages over conventional 
treatments. For example, one-particle dynamics of a gas is described by a Boltzmann equa- 
tion. One of the assumptions involved in the derivation of the Boltzmann equation is that 
the initial state of the system contains no correlations, which is clearly very restrictive 111 ]. 
However, in the general approach used here it is natural to break up any arbitrary initial 
state into a superposition of local density eigenstates, and that these may then be treated 
separately because of decoherence. The local density eigenstates typically have small or zero 
correlations. Hence, decoherence gives some justification for one of the rather restrictitve 
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assumptions of the Boltzmann equation. 

We have not estimated the degree of decoherence in the models considered here, although 
it could be estimated by looking more closely at the approximations involved in going from 
exact to approximate decoherence, described in Section 2. However, there may be a more 
rigorous (but more difficult) way of proving the results of this paper, which would allow the 
degree of decoherence to be estimated. This would be to prove a theorem similar to that 
proved by Omnes for phase space projectors [8|. For a system of N particles with phase 
space coordinates z = (p,q), Omnes considered (approximate) projection operators onto a 
region T (^> %) of phase space, defined by 

p r = j d N z |z)(z| (163) 

where the states |z) are some form of phase space localized states, such as coherent states. He 
showed that under certain reasonable conditions, the form of this projector is approximately 
preserved under unitary evolution, that is, 

e^ Ht P r e-^ m « P Tt (164) 

where T t is the original phase space region evolved along the classical phase space trajec- 
tories. It is easy to see that this ensures approximate decoherence of coarse-grained phase 
space histories and that the probabilities are peaked about classical phase space paths. The 
result is therefore very similar in spirit to the present paper. It seems very plausible that a 
similar result may be proved here for projections onto local densities. That is, we would like 
to construct a set of projectors onto the local densities, P ng h say, and then show that they 
are approximately mapped into P nt g t h t under unitary evolution, where n t , g t , h t are related 
to the initial values n, g, h by a closed set of evolution equations. Such a result is not simply 
obtained by a coarse-graining of the Omnes result, the issue being that n t ,gt,h t have to 
evolve according to a closed set of equations, which is not straightforward to accomplish in 
general. (The phase space coordinates p, q evolve according to the Hamilton equations, 
which is clearly a closed set of equations, but truncations or coarse-grainings of this set will 
generally not be closed). Moreover, the Omnes result breaks down when the underlying 
classical dynamics is chaotic. The corresponding hydrodynamic description, however, being 
coarser-grained, will generally not be chaotic and does not obviously break down, so this 
is a potential advantage of the hydrodynamic approach. These and related issues will be 
pursued in a future publication. 
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